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We present in this paper a detailed analysis of the flexoelectric instability of a planar nematic 
layer in the presence of an alternating electric field (frequency u), which leads to stripe patterns 
(flexodomains) in the plane of the layer. This equilibrium transition is governed by the free energy 
of the nematic which describes the elasticity with respects to the orientational degrees of freedom 
supplemented by an electric part. Surprisingly the limit lj — ► is highly singular. In distinct 
contrast to the dc-case, where the patterns are stationary and time-independent, they appear at 
finite, small oj periodically in time as sudden bursts. Flexodomains are in competition with the 
intensively studied electro-hydrodynamic instability in nematics, which presents a non-equilibrium 
dissipative transition. It will be demonstrated that ui is a very convenient control parameter to tune 
between flexodomains and convection patterns, which are clearly distinguished by the orientation 
of their stripes. 

PACS numbers: 61.30.Gd, 47.54.-r, 64.70.Md 



I. INTRODUCTION 



Nematic liquid crystals (nematics) are materials which 
prefer in some temperature range a uniaxial mean orien- 
tational ordering of their non-spherical molecules, while 
the positional order is fluid-like. The locally preferred 
axis is described by the director field n with n 2 = 1 0- 
3] . The basic equilibrium state, in the absence of external 
stresses, corresponds to a spatially uniform director con- 
figuration, where the orientational elastic free energy, F e i , 
with respect to n takes a minimum. Thus spatial varia- 
tions of n, which can be decomposed into splay, twist and 
bend distortions with the elastic constants kn, ^22, ^33, 
respectively, lead to an increase of F e i, i.e., to torques on 
n to restore equilibrium. In the presence of an electric 
field, E, and in non-equilibrium configurations electric 
torques on the director have to be considered as well. 
They originate from a nonzero electric polarization, P , 
which contains at first the standard dielectric contribu- 
tion Pdiei — £o( e — 1)E. The dielectric tensor, e, which 
depends on the local director orientation, is governed by 
the two dielectric permittivities eii and e± for E paral- 
lel and perpendicular to n, respectively; eo denotes the 
vacuum permittivity. Furthermore, in a rough analogy to 
the piezo-electric polarization of certain insulators due to 
mechanical strains, director distortions lead also to the 
so called flexo polarization, Pfi = ein(V ■ n) + e3(n ■ V)n, 
characterized by the two flexo coefficients e±, e^. If the 
electric field is strong enough the balance of electric and 
elastic torques may even require spatially periodic direc- 
tor variations where typically the flexoelectric torques 
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play a crucial role. The main goal of this paper is a 
detailed study of the resulting patterns, which are eas- 
ily detected by optical means due to the accompanying 
variations of the refraction index. 

Our analysis is restricted to the mostly used planar di- 
rector configuration, which is very convenient for exper- 
imental and theoretical studies. In this case a uniform 
basic state is achieved by sandwiching a thin nematic 
layer between two plates (parallel to the x — y plane), 
which may also serve as electrodes for the application 
of a transverse electric field E. By a suitable treatment 
of the plate surfaces the same in-plane director orienta- 
tion n — tiq — x is enforced at the confining plates, 
which then persists throughout the whole layer in the 
basic state. 

Before concentrating on the impact of the flexo torque 
we will briefly touch on the familiar consequences of the 
dielectric torque. It is destabilizing for a positive dielec- 
tric anisotropy, e a = en — ej_ > 0, in the planar case. Thus 
above a certain field amplitude the compensation of the 
dielectric torque by the stabilizing elastic one becomes 
impossible and the Freedericksz transition takes place: 
the director experiences a distortion in the x — z plane, 
maximal at the midplane of the nematic layer [l[ . The re- 
sulting director variation is in most cases uniform in the 
plane of the layer and varies only in the z-direction. Only 
in some exceptional cases, when kiil^w < 0.303, the ho- 
mogeneous Freedericksz transition for e a > is replaced 
by the so called splay-twist Freedericksz transition which 
leads to director distortions in the form of stripes parallel 
to no. We are only aware of one experimental realization 
[H, since almost exclusively 0.5 < ^22/^11 < 1 for ne- 
matics. Just for completeness, we mention that during 
the temporal evolution of the homogeneous Freedericksz 
state transient stripe patterns parallel to the x-axis have 
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also been described @. 

We now turn to the effect of a finite flexo polarization, 
which provides a robust mechanism (independent on the 
sign of e ) to generate stripe patterns (flexodomains) par- 
allel to no as well. The pattern forming instability of the 
basic state takes place when the flexo torque (oc |ei — e3 1 ) 
is sufficiently strong compared with the elastic and di- 
electric torques Q. Flexodomains are indeed observed 
in several nematic materials; also for e a < 0, where a 
Freedericksz transition is excluded 0-H[. In this paper 
we present a comprehensive analysis of the flexodomains 
in the presence of an electric field E = zEg cos(wi), both 
in the dc-case (circular frequency u) = 0) and the ac-case 
for small u. Unfortunately a previous analysis of the 
dc-case [ToL [Tl| suffers from a serious error. 

So far we have concentrated on equilibrium transitions 
driven by an electric field, which are governed by a min- 
imization principle of the free energy. In nematics, how- 
ever, we find often under the term electro-hydrodynamic 
convection (EC) a non- equilibrium, dissipative transition 
from the basic state towards a periodic arrangement of 
convection rolls. Since the resulting stripes run rather 
perpendicular to the initial director orientation no they 
are easily distinguishable from the flexodomains. The 
mechanism of EC, which has been first elucidated by Carr 
and Helfrich [l2|, [l3| in pioneering papers, has been com- 
prehensively discussed in a number of reviews in the last 
years (see, e.g., [13] and references therein). Some new 
aspects of the effect of flcxoelectricity on EC will be given 
in section ITVl of this paper. 

An interesting feature of our system is the possibility 
to switch between an equilibrium pattern-forming insta- 
bility (flexodomains) and the dissipative EC instability 
just by tuning the ac- frequency u. This scenario has 
indeed been observed in recent experiments on certain 
nematics Q: In the dc-case and at very small u one 
finds flexodomains, before EC takes over at increasing uj. 
Systematic theoretical analyzes, however, are practically 
missing so far. To provide those has been an important 
issue of our work. 

In section|TI]wc will briefly comment on the mathemat- 
ical background of our theoretical investigations. Sec- 
tion [HI] is devoted to the flexodomains. In section HVl we 
concentrate on EC at low frequencies in comparison with 
the transition to flexodomains. It will become obvious 
that in particular the limit uj — ¥ is far from trivial and 
requires a careful analysis. The paper concludes with 
some final remarks. 



II. MATHEMATICAL BACKGROUND 

Our starting point is a nematic layer (parallel to the 
x — y plane) of thickness d under the action of the applied 
potential U(t) — Eodcos(ujt) in the z-direction. In the 
quiescent basic state the director field is homogeneous 
throughout the layer (n = Uq = x) and flow is absent. 
Pattern forming instabilities are associated with a distor- 



tion, 5n, of no- They lead often to flow (v ^ 0) and to a 
perturbation, <fr, of the applied voltage U(t). Since both 
the planar director orientation at z = ±d/2 and U(t) are 
considered to be fixed, the conditions = and Sn = 
have to be fulfilled at the confining plates (z = ±d/2). 
This applies also to v under the realistic assumption of 
no-slip at the boundaries. 

The general theoretical framework to analyze electri- 
cally driven pattern forming instabilities in nematic liq- 
uid crystals is well established in terms of the standard 
nemato-hydrodynamic equations [U-HI . The electric field 
distribution is determined by the Maxwell equations in 
the electro-quasi-static approximation. As material pa- 
rameters we need the dielectric permittivities en, ej_, the 
electric conductivities <tii , o~± (for n parallel and per- 
pendicular to the electric field E), and the flexo coeffi- 
cients e\, e^. The director dynamics is driven by elastic 
and electric torques; in the presence of flow also viscous 
torques have to be taken into account. The velocity field 
is determined by the (generalized) Navier-Stokes equa- 
tion, where the viscous stress tensor depends on the ori- 
entation of n with respect to v and its gradients. For 
explicit calculations one needs the values of the elastic 
constants ka,i = 1,2,3 and furthermore the five inde- 
pendent viscosity coefficients oti,i = 1, . . . , 5 to quantify 
the stress tensor and the viscous torques. It is natural 
to introduce dimensionless material parameters of order 
one (labeled by primes). They appear then in the non- 
dimensionalized basic equations (see, e.g., [Uj]) and are 
usually defined as follows: 

ka = k' u ko , on = a^ao , 
(<r\\,o-±) = (c[|,Oj_)o- , 

(ei,e 3 ) = (e[,e' 3 )y/e k , (1) 

with 

fc = lCT 12 N , a = 1(T 3 Pa s , 

A s 

(T = 1(T 8 (fi m) \ e = 8.8542 x 1CT 12 - — . (2) 

V m 

For quantitative calculations in this paper we refer to 
the standard nematic MBBA [jq . which has been used 
in many experimental investigations in the past. The ma- 
terial parameters of MBBA are well known; for instance, 
for the analysis of flexodomains we use the following di- 
mensionless elastic and electric constants (the primes are 
omitted): 

ku = 6.66 , k 22 = 4.2 , 

ei = -3.25 , e 3 = -4.59 , e a = -0.53 . (3) 

By using Eq. (fTJ) it is conventional to measure lengths 
in units of d/ir, time in units of f, and to introduce a 
dimensionless control parameter R instead of the voltage 
amplitude Uq, where 

R _ e E$d 2 _ e E7o _ _ a Q d 2 

kpTT 2 kpTX 2 ' kpTT 2 
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To describe the onset behavior of the pattern form- 
ing instabilities the nemato-hydrodynamic equations are 
linearized about the basic state. As a result we arrive 
at a linear system of coupled partial differential equa- 
tions in the variables x, y, z, t for the perturbations 
5n = (0, n y , n z ), v, and </>. Since the lateral extensions of 
the nematic layer are much larger than the layer thick- 
ness d, periodic boundary conditions in the layer plane 
are appropriate. By switching accordingly with respect 
to the planar coordinates x = (x, y) to Fourier space 
q = (q,p), one gets the linear equations in the form to 
be found in the Appendix of 15*]; a simplified version to 
analyze the flexodomains will be presented explicitly in 
section IIIII 

As already explained, the perturbations Sn, v, <f> are 
assumed to vanish at z = ±<i/2; this is guaranteed by ex- 
panding these fields in terms of a complete set of Galerkin 
trial functions which vanish at z — ±n/2 in dimensionless 
units. For instance the director component n z (q,z,t) is 
represented as: 

M 

n z (q,z,t) = J~] n z (q,m,t)S m (z) , (5) 

with S m (z) = sin[m(z + 7r/2)]. We have tested that using 
a truncation parameter M = 4 yields already data to an 
accuracy of better than 0.1%; in many cases even M — 
2 is sufficient. Introducing the symbolic vector V(q,t) 
for the expansion coefficients n z (q,m,t), m — 1, . . . , M 
[Eq. ((SJ] and the corresponding ones for n y , v, </> we arrive 
after projection onto the trial functions at a linear system 
of coupled ordinary differential equations of the following 
form: 

C(q, t)^-V(q, t) = C{R, q, t)V(q, t) . (6) 

The matrices C and C are periodic with the ac-voltage 
period T = 2ir/ui. The general solutions of Eq. ([5]) have 
the Floquet representation V(t) — exp(at)Vo(t) with 
Vo(t + T) = Vo(t), where a defines the Floquet expo- 
nent. We are interested in time-periodic solutions V{t). 
Here a has to be purely imaginary of the form a = iuik/l 
with co-prime integers I > k. As a result the period of 
V(t) is given as l-T if k ^ and as T if k = 0. For a given 
q periodic solutions exist only for a discrete set of control 
parameters R — Ro(q) < Ri(q) < ^2(9) etc. Minimiz- 
ing Ro(q) with respect to q yields the critical wavevector 
q c and the critical control parameter R c = Ro(q c ) at 
which the quiescent basic state becomes unstable. In the 
context of the present paper the destabilizing modes at 
onset had always the period T, i.e., they are character- 
ized by a = 0. 

We have used two methods to calculate the periodic 
solutions of Eq. One option is to expand Vo{t) 

into a (truncated) Fourier series in terms of exp(inwt) 
with I n I < N. Then Eq. ^ transforms into an al- 
gebraic linear eigenvalue problem for the Fourier coef- 
ficients, from which we obtain Ro(q) (for more details 



see [15j). As will be documented below the time varia- 
tions in V(t) become increasingly sharper with decreas- 
ing w; consequently many Fourier modes up to N = 60 
had to be eventually kept. Thus, in an alternative, less 
time consuming approach, we construct numerically the 
matrix solution V{t) of Eq. © for the initial condition 
■P(O) = I where J denotes the unit matrix (see, e.g., [l7j 
and for a recent application [18]). We have to calculate 
the eigenvalues fix, fi2, . . . of the "monodromy matrix" 
M = V(t = T) where > |^a| etc- A periodic solution 
of Eq. © with period T exists when fix = 1. The small- 
est R to fulfill this condition yields again R = Ro{q)- 

As demonstrated in [15|], the general nemato- 
hydrodynamic equations as used in this paper, are in- 
variant against a reflection z — ¥ — z at the midplane 
combined with a translation in time by half a period 
T/2. Thus the solution manifold of Eq. (J6j) naturally 
splits into separate classes with different parity p — ±1, 
for which the notions "conductive" (p = 1) and "dielec- 
tric" (p = —1), respectively, have been introduced. For 
instance with respect to the director component n z one 
finds n z (— z,t + T/2) = pn z (z 7 t) (for more details, see 
[15|). This means that for the dielectric symmetry the 
time average of n z vanishes, while it is finite for the con- 
ductive symmetry. 

III. FLEXODOMAINS 

In the following we will investigate the bifurcation to 
flexodomains which leads to stripe patterns with the 
wavevector q — (0,p). It is easy to see that U(t) is 
not modified, i.e., = 0. Inspection of the full nemato- 
hydrodynamic equations (for instance in [HI]), shows 
that time variations of the director at nonzero oj lead 
in principle to a "back flow" which acts back onto the 
director in the form of viscous torques. They lead to 
small corrections to the dielectric and flexo torques of 
the order Oia^jrfe), with the Miesowicz coefficient r\i — 
(a 3 + a A + a 6 )/2 > 0. Since |a 3 /^ 2 | 2 = O(10~ 3 10~ 4 ) 
for MBBA and similar nematics, the viscous torques are 
safely neglected in this paper, which also facilitates the 
quasi-analytical approaches in section lTlI Bl below. More- 
over, this approximation has been validated by full nu- 
merical studies of the basic equations. 

It is convenient to introduce instead of the elastic con- 
stants fcn, fc 2 2 their average value, k av , and their relative 
deviation, Sk, from k av as follows: 

fcn = k av (l + 8k) , k 22 = k av (l - Sk) , (7) 

where obviously \Sk\ < 1. In contrast to the rod-like 
nematics like MBBA where £^2 < fcn, i-e., Sk > [see 
Eq. ||3J)], discotic nematics are characterized by fc 2 2 > 
fcn (Sk < 0) [lij]. Thus our analysis will cover negative 
Sk for completeness as well. Neglecting the back flow 
effects (v = 0) the director dynamics in flexodomains is 
only determined by dielectric and flexo torques. We start 
from the linear perturbations 5n(y, z, t) of the basic state 
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in position space and switch to Fourier space using the 
(real) ansatz: 



n y (y,z,t) = ain(py)n y (z,t) 
ra, (y, z, t) = cos(py)n z (z, t) 



(8) 



As a consequence the general equations © reduce to 
the following linear system of coupled partial differential 
equations (PDE's) in the variables z and t for the Fourier 
components n y and n z : 

d t n y = - [p 2 (l + Sk) - (1 - Sk)d zz ] fiy 

+ p [sgn(ei — ez)u cos(wi) — 26kd z ] n z , (9a) 
d t n z = - [p 2 (l - 5k) - fiu 2 cos 2 (ut) - (1 + Sk)d zz ] n z 

+ p [sgn(ei — es)u cos(wi) + 2Skd z ] n y . (9b) 

Instead of the time scale f and the main control parame- 
ter R [Eq. (dJ)] we have used in Eq. © the director relax- 
ation time Td and the dimensionless voltage amplitude it, 
which are defined as follows: 



1 e a R _ 1 £q e Q E 2 d 2 ^ 



with 



(ei - e 3 ) 



{ll/k av )f 
., and 7i = a 3 - a 2 ■ (10) 



Note that we will often refer in the following to the 
parameter fi as a reduced measure of the dielectric 
anisotropy e a . Furthermore it should be realized that 
in Eqs. ((SJ, (flO|) only the difference of the flexo coeffi- 
cients (ei — 63) comes into play. It is sufficient to confine 
oneself to the case (ei — 63) > 0, since the solutions of 
Eqs. ([9} for (ei — e 3 ) — > — {e\ — e 3 ) can be recovered by 
the transformation {n y (z), n z (z)} — > {n y (z), n z (— z)}. 



A. Flexodomains driven by a dc-voltage 

This subsection is devoted to the analysis of the flexo- 
domains in the dc-case (ui = 0) where Eqs. © are ex- 
act due to v = 0. Because all coefficients in Eqs. © 
are constant, the familiar separation ansatz fify tZ \(z,t) — 
e crt 7i(y ;Z ) (z) can be used to get rid of the time dependence. 
Thus the time derivatives d t in Eqs. © are replaced by 
a and one arrives at an autonomous linear system of or- 
dinary differential equations (ODE's) for the functions 
hr y<z \(z). As they have to vanish at z = ±7r/2, the eigen- 
values a will belong to a discrete set. The condition 
that the maximal eigenvalue a = ao(u,p) (growth rate) 
vanishes, determines the non-dimensional neutral curve 
uq(p). The minimum of uq(j>) at p = p c yields the criti- 
cal voltage u c — uq{p c ) 1 where (u c ,p c ) depend on 8k and 
/i. Alternatively uq(p) is given as the smallest u value, 
where Eqs. ^ have time-independent solutions, which 
vanish at z = ±7r/2. By the way, inspection of Eqs. (|9]) 
shows that uq (p) is even in p since it depends only on p 2 . 

Let us first concentrate on the case of a destabilizing 
dielectric torque with positive e a where /x oc e a > [see 



Eq. (TTOl) ]. For p = Eqs. ([9]) can be easily solved by 
choosing h y = and h z oc sin(z + tt/2). As a result we 
obtain 



u c (p = 0) 2 =u% = (l + Sk)/fi 



(11) 



This solution obviously describes the homogeneous splay 
Freedericksz distortion of the director. In fact, we recover 
from Up with the help of Eq. (p~0|) the familiar critical 
Freedericksz voltage Uq — Up = iTyJk Q k av (\ + 5k)/(e a € ) 
in physical units (which is not influenced by the flexo ef- 
fect). Clearly the condition uo(p c ) — u c < up is neces- 
sary for the prevalence of flexodomains with wavenumber 
p c ^ against the homogeneous Freedericksz distortion. 
It will turn out that they exist only for u less than an 
upper limit /-i max (5k), at which their critical wavenumber 
p c approaches zero. 

In the case of negative e a (u < 0), on the other hand, 
the dielectric torque is stabilizing. It overcomes eventu- 
ally the destabilizing flexo torques when fi approaches a 
lower limit fi m in(Sk) < from above, where p c diverges. 
Thus the director remains undistorted in the basic planar 
state for u < ^ m i n (Sk). 

In general, the z-dependence of the functions n y , z (z), 
which have to fulfill the ODE's introduced above, is cap- 
tured by an ansatz cx e Az . In our case we obtain four 
different values A = ±Ai, ±«A2, where only real A^ are 
compatible with the existence of flexodomains. The gen- 
eral solution, which consists of a linear combination of the 
four exponentials e ±AlZ , e ±lX2Z has to fulfill the boundary 
conditions of vanishing n y , n z at z = ±7r/2. As shown in 
the Appendix, one arrive thus at the following implicit 
equation for the neutral curve uq(p): 

A\ sinh(Ai7r) sin(A27r) 

+A2A1A2 [1 - cosh(AiTr) cos(A 2 tt)] = . (12) 

For the explicit expressions of the Aj, Ai, i = 1, 2, which 
depend on p, u, 8k, /i, we refer to the Appendix as well. 
Equation (|12p represents the neutral curve uo(p) in im- 
plicit form. Minimization of uq (p) with respect to p gives 
the critical wavenumber p c and the corresponding critical 
voltage u c = uq (p c ) of the flexodomains. 

In Fig. QJa) we show representative examples for u c 
as a function of /1 for Sk = and for Sk — ±0.3 cal- 
culated with the help of Eq. (TT^I) . The corresponding 
critical wavenumbers p c are shown in Fig. D^b). As dis- 
cussed before the Freedericksz state with p c = and 
« c (0) = up is smoothly approached when (Sk). 
Decreasing /1 from [i max on is associated with a mono- 
tonic increase of both u c and p c until they diverge at 
M = ^min{8k) < 0. It is obvious that the knowledge of 
the limit curves (J, m i n (Sk) and ^imax{8k) plays an impor- 
tant role to identify the regime of flexodomains in depen- 
dence on the parameters u and Sk. Thus we show these 
limit curves in Fig. [2]in the (u, Sk) plane, where \Sk\ < 1 
according to Eq. ([7]). On a first look, it is surprising that 
^max(Sk) diverges at Sk w 0.53. However, as discussed 
in detail in the following subsection, this divergence is 
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FIG. 1: (Color online) Critical voltage u c (a) and critical 
wavenumber p c (b) as a function of [i for Sk — (^ m in = — 1, 
Umax = 1), for Sk = 0.3 (pmin = —0.769, fj, ma x = 3.072) and 
for 8k = —0.3 (/^ m in = —1-427, ^ ma;c = 0.413). 



closely related to the existence of the (spatially periodic) 
splay- twist Freedericksz transition for 8k > 0.53 in the 
absence of the fiexo torque (ei — = 0) [4]. 

Before we turn, however, to further discussions of 
Eq. (|12p. we present at first a very useful approximate 
analysis of flexodomains for small \8k\. Instead of using 
Eq. (|12[) directly it is more transparent to start again 
from Eqs. @, by introducing the "one- mode" approxi- 
mation nt y>z -) oc sin(z + 7t/2). One obtains then immedi- 
ately at a quadratic equation for the eigenvalues a, the 
largest one determines the growth rate ao. Note that in 
this approximation the terms oc 2p8kd z rb( yz \ in Eqs. Q 
do not contribute. Standard perturbation analysis shows 
that they would produce corrections of the order 8k 2 to 
Co. The condition ao = leads within the one- mode ap- 
proximation to the following expression for the neutral 
curve: 



u I(p) 



(p 2 



1 



5k 2 (p 2 



1) 



p 2 + fi[p 2 + 1 + 8k(p 2 - 1)] 
with its minimum at p 2 = p 2 , where 



(13) 



The explicit expression for u 2 = u^(p c ) obtained from 
Eqs. (|T3|) . (fT4]l is quite lengthy and will not be shown. 
According to the general remarks above, the /i interval, 
where flexodomains exist, can in general be read off from 
p 2 given in Eq. (|14[) : the zero of the numerator deter- 
mines the upper limit p, ma x and the zero of the denomi- 
nator the lower limit fi m in ■ Thus we obtain the following 
approximate expressions valid for small 8k: 



/./,„ 



1 



l + 8k ' 



ASk) = 



l + 8k 
(1 - Skf 



(15) 



The one-mode approximation becomes exact in the spe- 
cial case 8k — (one constant approximation, fcn = ^22), 
where it yields a rigorous solution of Eqs. ©■ We recover 
in this case the results of fal: 



Pc 



1-M 
1 + M 



u 2 = 



(1 + M) 2 



(16) 



where /U TO i n (0) = -1 and p ma x(0) = 1- In general 
the exact neutral curve from Eq. (TT2j) is described to 
an accuracy of better than 0.5% by Eq. (fTS"]) for small 
\Sk\ < 0.2. We will demonstrate below, that /ij, in (ft) 
given in Eq. (|15[) even coincides with the exact curve 
Hmin(Sk) shown in Fig. ([2]) for all \Sk\ < 1. 

3 

2 

1 



=L 
-1 

-2 
-3 



homogeneous 
Freedericksz 



flexodomains 



basic state 




0.5 



FIG. 2: (Color online) Upper and lower limit curves, fi m ax{8k) 
(dashed, red) and fj,min{Sk) (solid, black), respectively, in the 
(Sk~n) plane. Hmax{Sk) diverges for 8k « 0.53, while fi m in — > 
— 00 for Sk — > —1. 

Finally we would like to stress that our analysis of 
flexodomains is at variance with recent investigations 
[l(J[ll| on the same subject. It will be explained in more 
detail in the Appendix that this work suffers from a ba- 
sic mathematical error. Thus for instance the prediction 
of a "singular" behavior of p c and u c both in the cases 
W fc 22 = 3 (Sk = 1/2) and k n /k 2 2 = 1/3 (8k = -1/2) 
does not hold. 



Analytical treatment of fj, ma x(Sk) and ti m in(8k) 



2 (-1 + 8k 2 )fi + V(l + Sk)[l + 8k(l + 4//)] 

Pc = Ti , n.NM , — , f ,-M ■ ( 14 ) 



(l + 8k)[l + ti(l + Sk)} 



So far we have given the exact description of flexo- 
domains through Eq. (|12p . In addition we have demon- 
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strated the usefulness of the analytical one-mode approxi- 
mation (I13p of the neutral curve u Q (p) at small 5k. In this 
section we will derive analytical expressions for the lim- 
iting curves fi m in(5k) and p m ax(5k) in the whole range 
-l<5fc<l. 

Let us start with the discussion of p m ax{5k) for the 
case /j, > 0, where we have competition between the 
flexodomains and the homogeneous Freedericksz state. In 
Fig.[3]we show a typical neutral curve uo(p) for 5k = 0.2 
and different /i obtained from Eq. ([12")) . At p = the 
function u (p) has an extremum with Uq(0) = up, see 
Eq. (jlll) . For /j, < [i ma x(5k) this point corresponds to a 
maximum where d 2 u^(p = 0) < 0; here and in the fol- 
lowing the notation 9™ for the derivatives (d n /dp n ) has 
been used. The minimum of Uo(p) at finite p = p c where 
dpUo(p — p c ) > and uo(p c ) — u c < up describes the 
flexodomains with wavenumber p = p c . 
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FIG. 3: (Color online) The neutral curve «o(p) normalized to 
the Freedericksz threshold uf as a function of p for 5k = 
0.2 and different p with pi = 1.4, p2 = 1.6, pa = 1.8, fJ.4 = 2.0. 

With increasing p both the critical wavenumber p c 
and the difference (up ~ u c ) decrease in Fig. [3] At 
A 4 = ^max{5k) = 1.965 the minimum and the maxi- 
mum of uo(p) merge at p = 0. Thus the equations 
u (p = 0) = up and d 2 uo(p = 0) = arc fulfilled, which 
can be solved with respect to \x by expanding Eq. (|12|) 
up to order 0(p 2 ). The resulting analytical solution is 
given as \i = p max (5k) with 



: (5k) 



l + 5k 



1 - 25k + 5k 2 (32/ t: 2 - 3) 



(17) 



It is convenient to introduce also the function 5kp(n) as 
the inverse of ~p max (5k), which is given as: 



5k F (n) 



l + 2p- y/{4p-l) 2 - 128/Lt(/2 - 1)/tt 2 
2^(32/tt 2 - 3) 



(18) 

Thus 5kp(pL) marks at fixed p the transition from the 
homogeneous Freedericksz state for 5k < 5kp(p,) to the 
flexodomains in the interval 5kp(p) < 5k < 1 (see also 
Fig. [3]). The special case 5kp(p = 1) = is consistent 
with p m ax(0) = 1 [see Eqs. (E 



From our reasoning it seems obvious that the ana- 
lytical expression ~p max {5k) should reproduce the curve 
^max{5k) shown in Fig. [2] As demonstrated in Fig. 0] 
this is indeed the case for the interval —5ki ow < 5k < 1 
where 5ki ow w —0.566. On this 5k interval we will first 
concentrate. 




FIG. 4: (Color online) Upper limit curve u max (solid line, 
black) as a function 5k from Eq. (|12[) in comparison with 
~p max from Eq. fTf) (dashed line, red). 



One sees immediately that 
V-max^k) approaches zero (i.e., JI max 



approaches a critical value 5k c 
5k c is given as: 



the denominator of 
diverges) when 5k 
> from below, where 



1 - 2Jl -8/tt 2 
5k c = 32 V a _ 3 / = 0.5346 



(19) 



Consistently 5kp(fi) in Eq. (1181) approaches for /i — > oo 
the limit 5k c . As a result, flexodomains are exist for any 
/i > when 5k > 5k c , while they are restricted to the 
region p < ~p max {5k) in the case of 5k < 5k c . 

It is not a coincidence that 5k c is very near to the num- 
ber 5ksT ~ 0.53, which has been quoted in the literature 
more than two decades ago in a different context: Ac- 
cording to |4j the homogeneous Freedericksz transition, 
in the absence of flexo effects (e.\ — e% = 0), is replaced 
for 5k > SksT by the spatially periodic splay-twist (ST) 
Freedericksz transition. Similar to the flexodomains a 
state of finite n y , n z bifurcates from the basic state which 
is periodic in the y-direction with a critical wavenumber 
p c , to balance the dielectric and the elastic torques. In 
this case additional flexo torques should not play a crucial 
role and the striped director configuration should develop 
even at arbitrary small \ei — ea\, which corresponds to an 
arbitrary large p = e a k a v/{ei — e^) 2 . In fact, as shown in 
the Appendix, the limit ei — e3 = is covered by Eq. (TT2|) 
and one finds Sksr = 5k c . 

Our theoretical considerations are confirmed by the 
representative numerical results for the critical voltage 
u c and the critical wavenumber p c of flexodomains, which 
are shown in Fig. [S] as function of 5k < 1 for two differ- 
ent fi. In the case fi = 2 the homogeneous Freedericksz 
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state is replaced at 5k — 5kp(p, = 2) « 0.206 by the 
flexo domains; with increasing 5k the critical voltage u c 
monotonically decreases from u c — up on, while p c in- 
creases from p = on. For — > oo (i.e., in the absence 
of flexo effects) with 5k p — > 5k$T we get in the interval 
5ksT ^ 5k < 1 the splay-twist Freedericksz distortion as 
a special case of the flexo patterns. 




FIG. 5: (Color online) Relative critical voltage u"i/u 2 F and the 
critical wavenumber p c as function of 5k for p, = 2 (solid line 
for u^/up and dotted line for p c , black) and in the absence 
of flexo effect (ej — ez = 0) (dashed line for v%/u% and dot- 
dashed line for p c , red). For p, — 2 the Freedericksz state is 
approached at ~ 0.206, while this happens for ei — e-i — 
at <5fc = 5fcsT [data from Eq. ijH) ]. 



We will now return to the case — 1 < 5fc < <5/c; ou) s» 
—0.566 in Fig. [4] where the exact fi m ax(5k) from Eq. (fT2|) 
is slightly larger than ft m ax(5k) given in Eq. (IT71) . The 
reasoning, which led to the expression pL max , is rigorous 
as long as the p-dependence of neutral curves is of the 
type shown in Fig. [3] As demonstrated for instance in 
Fig. [5] for 5k = —0.7 this does not hold in the vicinity of 
5k = —1. The condition uq(p) = up, which determines 

t^rnax 

(5k), is now already realized at a finite p = p c , where 
uq(p) has a second local minimum of height uo(p c ) = up. 
Since a merging of the two extrema is not required, as 
assumed in the calculation of p(5k) from Eq. ([T7]) . we 
have now fi m ax(5k) > ~p max (5k) in agreement with Fig.@] 
The exact neutral curve for 5k < 5ki ow is only accessible 
by numerically solving Eq. (|T2"1) . The special value 5k = 
5ki ow , where \x max starts to deviate from p m ax{5k), is 
obviously determined by the condition d^uoij) = 0) = 
dpUo(p = 0) = at u = up and /i = Hmax(5k). With 
the use of Eq. (fT2|) one obtains from these conditions the 
following closed analytical expression for 5ki ow : 



5k[ ow — 



tt 2 (96 - 5tt 2 - 8Vtt 4 - 54tt 2 + 468) 



-137r 4 + 8327r 2 -6912 



= -0.5666 



1/2 



(20) 



which agrees perfectly with the results presented in 
Fig. H 

Finally we consider the case p, oc e a < 0, where the 
dielectric torque is stabilizing and where the competing 
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FIG. 6: (Color online) The neutral curve «o(p) as a function 
of p for 5k — —0.7 and different \x where ui = 0.11, = 
0.12, /Lt 3 = 0.13, /i 4 = 0.14. 



Freedericksz transition is absent. As already mentioned 
we found fJ, m in{5k) = ^ in (<5fc) given in Eq. (fl~5j) . This 
can be explained by the fact, that exactly at (J, = p a min 
the purely imaginary roots ±?A2 in Eq. (|12p become real 
(see the Appendix). As a consequence the trigonometric 
functions sin(A27r) and cos(A27r) change into hyperbolic 
ones; then Eq. (fT2]l allows only for the trivial solution 
u = 0, p = 0. 



B. Flexodomains driven by an ac- voltage 

While the onset of flexodomains in the presence of 
a dc-voltage is well described by the closed expression 
Eq. (fT2"|) , the case of an ac- voltage requires in general nu- 
merical methods along the lines described in Sec. [TTJ We 
consider only low frequencies u where ujt^ < 20 which 
corresponds to frequencies / = uj/(2tt) up to 20 Hz for 
a cell of thickness d = 10 jum filled with MBBA. Note 
that we will sometimes use in the ac-case instead of the 
voltage amplitude, Uq — Egd, the effective (rms) voltage 
f/o/%/2, which derives from the time average of U 2 (t). 

As a warm-up we first consider the Freedericksz tran- 
sition ((!«£„> 0) at finite uj, where n y = and p = 
in Eqs. ([9]). In analogy to section UlI Al we obtain the 
solution n z (z, t) — n z (t) sin(z + 7r/2) with: 

n z {t) = n z (0)exp|-^ dt'[l + 5k - fiu 2 cos 2 (ut')} j (21) 

The solution h z (t) is periodic with period T = 27t/u 
when the integral in Eq. (I2ip vanishes for the upper limit 
t = T. Thus the Freedericksz transition voltage, «j? r (w), 
for uj ^ is obtained as: 



u Fr (uj)/V2 =(l + 5k)/fx 



(22) 



Obviously it is the (non-dimensionalized) effective (rms) 
critical voltage Up r (uS)/ \/j2), which continuously ap- 
proaches in the limit uj — >■ the Freedericksz transition 
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voltage up (see Eq. (|TT1) in the dc-case. At onset for 
u = up r the solution Eq. (|2ip. normalized to its maxi- 
mal value, has the following form 

f , r d - sin(2wi) 1 , , 

n z {t) = exp | -(1 + 5k) ^ | . (23) 

At small oj Eq. (|23|) represents sudden burstdike director 
distortions similar to those shown in Fig. |H1 The peaks 
of the h z (t) are located at tut — 7r/4, 57r/4, etc. and then- 
width is proportional to y/Zj. 

Now we turn to the discussion of flexodomains where 
n y 7^ and p / 0. The threshold u c and the criti- 
cal wavenumber p c for finite frequency uj are calculated 
as explained in Sec. [TT] As a rule we get two types of 
solutions, one with the conductive, the other with the 
dielectric temporal symmetry. In the following we use 
in general the material parameters of MBBA given in 
Eq. ©, where k av = 5.43, 5k = 0.2265 [see Eq. 0] and 
ei — e3 = 1.34. The dielectric torque is stabilizing, since 
e a = —0.53. According to our investigations the existence 
of flexodomains in the dc-case seems to be in general a 
necessary prerequisite for their existence in the ac-case. 
In fact, at u = flexodomains solutions do not exist 
for standard MBBA, since according to Eq. (fTUf we have 
/j = —1.6027 < Hmin{5k) = —0.815; they have also not 
been observed in experiments. Thus we follow an idea in 
[2fj| and replace the difference (ei — ea) of the flexo coef- 
ficients in Eq. © by the product £_ • (ei — 63) to study 
specifically the impact of the flexo torque. Flexodomains 
then exist for > 1.402. 

Representative examples for the dependence of the 
critical voltage u c and the critical wavenumber p c on the 
"flexo strength" £_ and on the ac-frequency w (in units 
of t~[ ) are shown in Fig. [71 For £_ = 3, corresponding 
to /1 = —0.178, the conductive symmetry prevails in the 
whole frequency range, while for £_ = 2, corresponding 
to fi = —0.401, a switch over from the conductive sym- 
metry to the dielectric symmetry takes place at larger 
LOTd- For both symmetry types the u c and p c curves rise 
monotonously as function of uiTd at fixed On the 
other hand, for the same symmetry type u c and p c are 
seen to decrease with increasing i.e., with growing 
magnitude of the flexo torque in analogy to the dc-case. 

To study the dependence of the critical voltage u c (ui) 
and the critical wavenumber p c (w) on the elastic con- 
stants, still MBBA parameters have been used, except 
a change of 5k. Furthermore we have chosen £_ = 2 
to favor the flexodomains against the homogeneous basic 
state. The resulting data are shown in Fig. [5] As before 
u c and p c increase monotonically with lot^- We are, how- 
ever, unable to offer a simple explanation for the change 
of symmetry with 5k: While the dielectric symmetry of 
the flexodomains solutions prevails for 5k = —0.3, 0, we 
find the conductive symmetry for 5k — 0.3. The switch 
over between the symmetries happens at 5k sa 0.1. 

In Fig. [S] we show an example of the director dynamics 
in the midplane (z = 0) as a function of time at the low 
frequency UTd — 0.05. Remarkable are the sudden burst- 
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FIG. 7: (Color online) Critical voltage u c (u) in dimensionless 
units (a) and critical wavenumber (b) of the flexodomains 
as function of WTjj for MBBA parameters with £_ = 2 and 
£_ = 3. The symmetry of the solution is indicated. 



like director distortions, i.e., patterns would appear in 
the experiments only for very short time intervals. The 
conductive time symmetry for 5k — 0.3 is reflected in the 
finite time average of n z (t) ; it vanishes for n y (t) . On the 
other hand, for 5k < 0.1 the behavior of n z and n y with 
the dielectric time symmetry would be just opposite. 

Before we have concentrated on the MBBA-like ma- 
terials with fi < 0. In the case of fj, > we are faced 
with the competition of flexodomains with the Freeder- 
icksz transition as discussed in Sec. lIII Al for dc-case. The 
threshold voltage for the flexodomains increases mono- 
tonically with u) similar to the case of negative \x (see 
Figs. [JJ Since on the other hand u Ft {uj) in Eq. (|2"2j) 
does not depend on w, the Freedericksz transition will re- 
place the flexodomains in any case above a certain finite 
lo. The frequency range where the flexodomains prevail 
will shrink with increasing fi. 
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FIG. 8: (Color online) Critical voltage u c (uj) in dimensionless 
units (a) and critical wavenumber (b) of the flexodomains as 
function of UTd for different 5k and £- = 2. 



The limit to — > 

The limit lj — > is obviously not smooth: In the dc- 
case both n z and n v have finite amplitudes at onset, while 
in the ac-case the time-average of one of the director com- 
ponents vanishes depending on the time symmetry. An 
interesting question is, how this discontinuity is reflected 
in the critical voltage. As a rule of the thumb it is typ- 
ically assumed, that the effective (rms) voltages in the 
limit lj — > do not differ too much from critical voltages 
in the dc-case. For the Freedericksz transition that differ- 
ence is even zero [see Eq. (|22p ] and for electroconvection 
it is of the order of 10% (see the Appendix of [H|). In 
the case of flexodomains the corresponding differences 
are considerably larger in most cases. 

As a representative example we discuss the w — > limit 
of the u c , p c curves in Fig. [7] in comparison with the 
dc-case from Eq. (fl2l . Here MBBA material parameters 
have been used except that we allow for a modification of 
ei — e3 via a factor £_ introduced above. At low frequen- 
cies the threshold curves both for £_ = 2 and £_ = 3 have 
conductive symmetry. For £_ = 2 (fi = —0.401) we find 
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FIG. 9: (Color online) The director fields as a function of time 
at flexodomains threshold for a periodic excitation. uiTd = 
0.05, 5k = 0.3 and £_ = 2 (conductive symmetry). 



for instance u e //(0) = u c (lj — > 0)/y/2 = 4.52 compared 
to the dc-value u c (0) = 3.49; for = 3 (fi = -0.178) 
the difference between u e ff(0) = 2.83 and u c (0) = 2.42 
is smaller. The discontinuity of u e ff(cj) at lj = is asso- 
ciated with corresponding discontinuities of p c : we find 
p c (u! —> 0) = 1.66 compared to p c (0) = 1.5 for £_ = 2 
and p c (u —> 0) = 1.25 compared to p c (0) = 1.17 when 
£_ = 3. For the parameters chosen in Fig. [8] we find 
comparable discontinuities as in Fig. [7] for = 2, where 
the symmetry of the director configuration (dielectric for 
8k = —0.3, and conductive for 6k = 0.3) does not seem 
to play a significant role. 

The discontinuities in the limit to — > can be demon- 
strated most clearly in the special case 6k = 0. Then 
the z-dependence of the director fields in Eqs. ([9]) is 
rigorously described by the ansatz {n y (z,t),n z (z,t)} = 
{n y (t), h z (t)} sin(z + 7r/2). The transformation of the 
time variable t — > t/u) leads to: 



ujdthy = — (p 2 + l)n y + sgn(ei — e^)pu cos(t)h z 
ojd t n z = — [p 2 + 1 — /iM 2 cos 2 (i)]n z 
+sgn(ei — es)pu cos(t)n y . 



(24) 



From a mathematical point of view Eqs. (|24p present an 
interesting dynamical system. It is kind of "singular" 
since the prefactor of the time derivatives is proportional 
to lj and thus small in the limit u) — > 0. For the special 
case fx = (e a = 0) the general solution of Eqs. (|2"4")l is 
given as: 



hy (t) = cxp 
n z (t) = exp 
with <p{t) 



p 2 + 1 



cie 0(*) _ C2e -^( 



p 2 + l 



"}• 



sgn(ei - e 3 )pusin(i) 



(25) 



where c\, ci denote integration constants. Inspection of 
Eqs. (|2"5")) shows immediately, that one is unable to find 
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values of c\, ci to allow for periodic solutions. This result 
is a further hint, that the limit U) — > is not trivial, since 
in the dc-case according to Eqs. (fl6|) flexodomains do 
exist at /i cx e a = 0. 

The system Eqs. (|24|) for /i ^ has been in general 
investigated by the methods presented in Sec. El To 
study in particular the limit lo —¥ we have employed 
the approximation scheme presented in the Appendix of 
[2lj . which is closely related to the familiar WKB ap- 
proximation [22| in quantum mechanics. We use the 
ansatz {n y (t),h z (t)} = {h y (t), n z (t)} exp[S , (i)/w] and 
neglect the derivatives of h v {t), n z (t). The resulting cou- 
pled homogeneous equations for n y (t), n z (t) will have a 
nontrivial solution when the corresponding determinant 
vanishes. In this way the following equation for So(t) is 
obtained: 



{dtS f + gi (t)d t So + g (t) =0 



(26) 



where 



gi{t) = 2(p 2 + 1) - fiu 2 cos 2 (i) , 

g (t) = ( P 2 + l) 2 - [ P 2 + MP 2 + I)]" 2 cos 2 (t) . (27) 



From Eq. (|26p we obtain: 



2d t S 



-91 



9l 



4 5o 



(28) 



In analogy to the discussion of the Freedericksz transition 
[see Eq. (|2ip] the solution is bounded when the average 
over one period (2ir) of the right hand side of Eq. ([28]) 
vanishes. The required time integration can be performed 
analytically and leads to the following implicit equation 
for the neutral curve uq(p) = Mo(p, oj — > 0) 

n[2(p 2 + 1) - /i M 2 /2] \f,\ = 

2^V^ + (V + 11 2 Uq) arcsin (sj A /f^ u ^ ) 

This transcendental equation is solved numerically, which 
yields within the WKB approximation the critical volt- 
age and the critical wavenumber p^ of the flexo- 
domains in the limit lo — > 0. The data agree perfectly 
with the exact solution of Eq. (pM)) on the basis of the 
monodromy matrix method in the limit lo — > 0. In 
Fig. [TUTa 1 ) the ratio of dimensionless effective (rms) criti- 
cal voltage (it^/v2) and the dc- value is plotted as func- 
tion of /i. While this ratio is near one for — 1/2 < /j, < 1 
it increases strongly when further decreasing fi towards 
fi = /i m i n (5k = 0) = — 1. The corresponding ratio of 
the critical wavenumber p^ and the dc- values is shown 
in Fig. [TUTb). In comparison to the critical voltages one 
finds even larger deviations of this ratio from one in par- 
ticular at larger positive fi where p c itself becomes small. 

Finally we would like to point out, that the transi- 
tion to the dc-limit and the interesting time evolution 
of director perturbations characterized by bursts depend 
also on the wave-form of the exciting ac-voltage. To 
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FIG. 10: The critical effective voltage for flexodomains 
«c /V2 (a) and the critical wavenumber pY (b) from Eq. (129p 
normalized to the dc-values. 



in the limit uj 



demonstrate this we have considered a square- wave exci- 
tation which consists of an alternating sequence of con- 
stant voltages ±u on time intervals of the length T/2, 
i.e., U(t) — [/osgn[cos(wi)]. Then Eqs. (O can be solved 
quasi-analytically by joining continuously the analytical 
solutions on the parts with constant voltage, which leads 
to an implicit transcendental equation ("boundary de- 
terminant") for the neutral curve uo(p). Since for small 
ojTd the values of u c and p c are practically determined 
by the long-time intervals of constant u they are practi- 
cally identical to the dc-values. Consequently, as demon- 
strated in Fig. 111! the time variations of n z , n y are rather 
smooth compared to the bursts observed with the har- 
monic ac-driving at low lo. 

In the following section, we will study the alterna- 
tive electroconvection (EC) solutions of the nemato- 
hydrodynamic equations, which compete with the flexo- 
domains. In general it will turn out, that only for small 
to the flexodomains have the chance to prevail. 



IV. ELECTROCONVECTION 

The salient elements of the positive feedback loop 
leading to the dissipative electro-hydrodynamic instabil- 
ity (EC), which are contained in the familiar nemato- 
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tive (rms) critical voltages in physical units everywhere 
in this section. 



FIG. 11: (Color online) The director fields as function of time 
at flexodomains threshold for a square- wave excitation. u)Td — 
0.05, Sk = and £_ = 2 (dielectric symmetry). 



hydrodynamic equations (see Se c. ITfl) . have been eluci- 
dated by Carr and Helfrich [l2|, Il3j : A necessary con- 
dition is that the nematic is "contaminated" by a small 
amount of mobile ions, which results in a finite, though 
very small electrical conductivity of the order of er = 
1CP 8 (fi m) -1 . Then spatial fluctuations of n in the 
presence of a nonzero E lead to an electric charge den- 
sity p e i = V • D + pfi. The first contribution to p e i is the 
standard one, which derives from the dielectric displace- 
ment D — eoeE. The latter "nonstandard" contribution, 
the flexo charge pji, is determined by the earlier intro- 
duced flexo polarization as follows: Pfi = V-P//. Via the 
Coulomb force, p e iE, in the (generalized) Navier-Stokes 
equation a material flow is driven, v, which exerts an 
additional viscous torque on the director. Under a favor- 
able constellation of the material parameters the viscous 
torque reinforces the initial director distortion leading to 
EC. For more details in particular a most recent review 
[l~4j might be useful, where systematically the sensitive 
influence of the sign of the dielectric anisotropy e Q , of 
the electric conductivity, er a , and of the basic director 
configuration on the patterns is discussed. 

Exploiting the calculational scheme briefly discussed 
in section [TT] we have precisely characterized the onset of 
EC. While in flexodomains the time scale is set by the 
director relaxation time oc d 2 [see Eq. (|10p], in electro- 
convection the thickness independent charge relaxation 
time T q = eoe± /crj_ plays an important role as well. Thus 
we observe in EC a thickness dependence of the critical 
properties, which cannot be absorbed by a suitable fre- 
quency rescaling. For definiteness we concentrate in the 
following on a nematic layer of thickness d — 10 pro. and 
use MBBA parameters (if not otherwise stated), where 
Td = 0.2 s and r q = 4.65 x 10 -3 s. The impact of the flexo 
charge pji on EC is determined by the sum (ei+ea) of the 
flexo coefficients. In analogy to the factor introduced 
in section IIIII to tune the strength of the flexo torque, 
we replace here in addition (ei + e 3 ) by £ + ■ (ei + 63) to 
allow for a selective modification of /?/;. To facilitate a 
direct comparison with experiments we will refer to effec- 



In Fig. [T2] we present representative critical EC-data 
for MBBA parameters [l6| as function of uir q for differ- 
ent £ = = £ + . In general U c and \q c \ [see Figs. [T27 a). 
(b)] increase monotonically with increasing u>. While the 
two curves are smooth for £ = £_ = = 2, where the 
dielectric symmetry prevails, we observe discontinuities 
for the two smaller £. Here the conductive branch at 
small u) is replaced by the dielectric one at a crossover 
frequency u> c where uj c T q s» 1.2. The obliqueness of the 
rolls, as shown in Fig. [T2T c). is measured by the angle a 
between q c and the preferred director orientation in the 
basic planar configuration (our x-axis) . In the absence of 
flexo polarization (£ = 0) we find a = for all ui. When 
= £ + = 1 the obliqueness of the rolls in the conduc- 
tive regime vanishes continuously at the "Lifshitz point" 
uJLT q ~ 0.35 and becomes finite again in the dielectric 
regime. At the largest £ — 2 the angle a ~ 40° is finite 
for all oj. Note that recently for other nematic material a 
"Lifshitz point" in the dielectric regime of EC has been 
identified as well [9(. 



In analogy to the flexodomains the limit u> — > for the 
critical voltage in EC is not smooth. However, the dc- 
critical voltage is much better approximated by the limit 
u) — > of the effective voltage U c {uj) than for the flexo- 
domains. This finding can be explained within a WKB 
approximation [2l| , where the relative corrections are in- 
deed small of the order of T q /rd (up to some g-dependent 
factors). EC shares, however, with the flexodomains the 
typical spiky time evolution of the director components 
at very low oj, as comparison of Fig. [5]with Fig. 1 131 shows. 
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FIG. 12: (Color online) Critical effective voltage U c of elec- 
troconvection (rms value in volt) (a), modulus of the critical 
wavevector |q c | (b), and the angle a of the critical wavevector 
with the a>axis (c) as a function of frequency oj (in units of 
the charge relaxation time r q ) for different magnitudes of the 
flexo coefficients £ = = (for more details, see text). 



Besides the fiexo coefficients a further crucial parame- 
ter is the anisotropy a a of the electric conductivity, since 
the Carr-Helfrich charge separation mechanism depends 
strongly on the magnitude and the sign of a a /a±. This 
dependence had not been investigated in detail for many 
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FIG. 13: (Color online) The director fields as a function of 
time at EC threshold for u)T q = 2.3 x 10" 3 (ujT d = 0.1). 



years in the past, since in the materials used at that time 
c a /oj_ would only vary between 0.3 and 0.7. Recently, 
however, materials have been found, where by decreas- 
ing the temperature a a decreases as well; it passes even 
zero and becomes negative. Thus it has now become 
important and attractive to study the a a dependence in 
more detail. This has been done partially in [lj|, where 
U c (uj) has been observed to move up in general with de- 
creasing a a - To study the details we have taken MBBA 
parameters as before except that cr a /a± was allowed to 
vary. Furthermore, to explore the competition with flexo- 
domains, which do not exist for standard MBBA (see 
Sec. IIII B|) the flexo coefficients have to be increased as 
well. We have centered on the case £ = = £ + = 2, 
which has been shown already in Fig. [12] for EC and in 
Fig.[7]for flexodomains. The EC threshold, which is char- 
acterized by the dielectric symmetry, increases strongly 
with decreasing a a /cr± for all oj as documented in Fig. 1141 
Only for (T a /a± = —0.5 and for very small frequencies up 
to ujT q » 0.1 (where ujt^ s» 3.6) flexodomains with con- 
ductive symmetry prevail. For the case = £+ = 1 
one finds the same behavior of U c (uj) when changing a a , 
while the flexodomains do not exist. 

To find material combinations, where flexodomains 
have a chance to dominate the pattern forming instability 
for standard and not too small values of a a , one has ob- 
viously to chose different values of £_ and £ + . One needs 
a relatively large £_ to increase the flexo torques, which 
are responsible for the existence of flexodomains. Con- 
versely, £+ should be fairly small to suppress the positive 
influence of the flexo charge on EC. Since systematic ex- 
perimental studies, also including different nematic mate- 
rials are missing so far, extended parameter studies seem 
to be futile at the moment. So we present only one ex- 
ample with £_ = 3, £ + = 0.25 and a a /a± = 0.2 where 
flexodomains (with conductive symmetry) indeed prevail 
EC at low frequencies as shown in Fig. [151 
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FIG. 14: (Color online) Critical voltage U c (rms value in volt) 
for EC as a function of frequency u> (in units of the charge 
relaxation time r q ) for five different values of a a /a± between 
0.5 and —0.5 and for f_ = = 2. The corresponding critical 
voltages for the flexodomains (conductive symmetry at small 
lu and dielectric symmetry at larger cu) are included as well. 
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FIG. 16: (Color online) Critical effective voltage U c of elec- 
troconvection (rms value in volt) as a function of frequency uj 
(in units of the charge relaxation time r q ) for MBBA material 
parameters (£+ = £_ = 1) except a a /<J± — —0.05. 



Finally, we address briefly the intriguing limiting be- 
havior of EC at low frequencies, when £_ is too small for 
the existence of flexodomains at u> = 0. As demonstrated 
in Fig.[T5]we predict in this situation a persistent switch- 
ing between the conductive and the dielectric solutions 
for <J a /v± = —0.05. 



V. CONCLUSIONS 
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FIG. 15: (Color online) Critical effective voltage U c (rms value 
in volt) as a function of frequency ui (in units of the charge 
relaxation time r q ) for flexodomains and EC patterns. Mate- 
rial parameters of MBBA except a a /cr± — 0.2, £_ = 3, and 
£+ = 0.25. 



In this paper we have investigated pattern forming in- 
stabilities in nematic liquid crystals originally in the basic 
planar configuration, which are driven out of equilibrium 
by an electric dc or ac voltage applied perpendicular to 
the layer. In the basic state the director field is homoge- 
neous over the nematic layer and a preferred direction is 
singled out in the layer plane. 

In particular we have concentrated on the impact of the 
flcxo effects, which generically exist, when the director 
varies in space. They results in on one hand to the fiexo 
torque which may lead at sufficiently large voltage ampli- 
tude to the flexodomains. These are spatially periodic in 
the layer plane and can be visualized in experiments as 
stripes parallel to the preferred planar direction. In this 
paper the bifurcation scenario to flexodomains is for the 
first time comprehensively analyzed for applied dc and 
ac voltages. We have also demonstrated that the well- 
established splay-twist Freedericksz transition in the ab- 
sence of flexo effects can be described in this framework. 
Particular focus was on the intriguing phenomena in the 
limit of vanishing ac-frequency w of the driving voltage. 
One finds for instance periodic-in-time, burst-like excur- 
sion of the director from the planar basic state. 

Furthermore we have investigated the competition be- 
tween the flexodomains and the familiar EC convective 
rolls. The latter are visualized as stripes which include a 
finite angle typically larger than 45° with the preferred 
axis. In general EC prevails, except at very low fre- 
quencies where the flexodomains may trigger the desta- 
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bilization of the homogeneous basic state. Like the flexo- 
domains the EC pattern appear also periodically in time 
as sudden bursts, while the system remains otherwise in 
the structureless basic planar state. 

We expect that our theoretical analysis will moti- 
vate further investigations of electrically driven pattern 
formation in nematics with special focus on the low- 
frequency regime. Though in particular the sudden, 
periodic-in-time "blow-up" of flexodomains and EC pat- 
terns as well a crossover between them have been clearly 
seen in recent experiments a quantitative compari- 
son with the theory is far from trivial. On one hand 
one is faced with considerable uncertainties in the ma- 
terial parameters even for a nematic like MBBA. More- 
over, one has to realize that the nematic layer is con- 
fined by metallically coated glass plates with a thin poly- 
meric alignment film on top, which themselves may have 
fairly complicated electric properties. Thus the whole 
system is represented as an equivalent circuit diagram 
(see, e.g., [H, where the total voltage applied to the 
cell is different from "theoretical" results like U c , which 
represent only the voltage drop over the nematic layer. 
The necessary corrections to U c are not easy to model 
and would also require systematic measurements on the 
empty cell. Finally the coexistence between flexodomains 
and EC pattern in the nonlinear regime when their linear 
thresholds are comparable has to be investigated. For in- 
stance such a scenario has indeed been observed in hybrid 
aligned nematic cell [25| in the presence of a dc-voltage. 
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Appendix A: Flexodomains driven by a dc-voltage 

In this Appendix we describe briefly the derivation of 
the implicit equation (fT2"j) which yields the neutral curve 
uo(p). We use the ansatz {n z (z),n y (z)} — {n z ,n y }e Xz 
in Eq. and set a — from the beginning. Thus we 
arrive at a linear homogeneous 2x2 system for n y n z 
which leads to the following secular equation for A: 



The eigenvector of the linear system is chosen as: 

r- f\\ - f\\\ Ji sgn(ei - e 3 )w - 2(5fcA 1 
K(A), % (A)} = |1, + gk) _ (1 _ 6k)x2 | .(A3) 

In line with the standard procedure we express the gen- 
eral solution as a linear combination of the eigenmodes 
[1 



[p 2 {l + Sk) 
-P 2 (sgn 2 (ei 



(l-Sk)X 2 } [p 2 (l 
- e 3 )u 2 - 45k 2 X 2 ) 



5k) 
. 



liu 2 - (1 + 8k)X 



3 — Xiz 



3 — iA 2 z 



n z {z) = aie"- 1 " + a 2 e + 036"""" + 04c 
n y (z) = aih y (\i)e XlZ + a 2 h y (-Xi)e~ XlZ 

+a 3 h y (i\ 2 )e lX2Z + a 4 h y (-i\ 2 )e- lX2Z , (A4) 

to fulfill the boundary conditions n z (±7r/2) = 
n y (±ir/2) = 0. A nontrivial solution exists, when the 
determinant of resulting set of four linear homogeneous 
equations for the a, vanishes. The resulting implicit 
equation for the neutral curve has been already given be- 
fore in Eq. (fl"2j) . It contains the \ (|A2[) and coefficients 
Ai, which are given as follows: 



Ax 



(f 2 1 x 2 ~f 2 x 2 )c 2 + (x 2 + x 2 YD 



with 



A 2 = 2hf 2 C 2 , 



h=X 2 -B, f 2 = X 2 + B, B=p 



(A5) 



1 + 5k 
I -8k ' 



C = 



25k 



VI - 5k 2 



D = 



sgn(ei - e 3 )u 
VI - 5k 2 



(A6) 



From Eq. (|A1|) we obtain the two eigenvalues A 2 = 
A 2 , —A 2 ., where A 2 , A 2 . > 0, which read as follows: 



\ 2 2 
X x =p - 



2(\ + 8k) 



sgn 2 (ei 



e 3 j 



I -5k 2 



S, X\ = -p 2 



25kfi 2 2 
p u 



l±u 



2(1 + 8k) 



2(1 + 8k) 



+ 5, 
1/2 
(A?) 



Maple and Mathematica have been very useful to per- 
form and to validate all our calculations. 

The Freedericksz transition in the absence of flexo ef- 
fects (ei — e 3 — 0) is formally covered by Eq. ([T2"j) as 
well. First of all the voltage u appears always in the 
combination (nu 2 ) where according to Eq. (JTOJ) the fac- 
tors (ei — 63) cancel. Furthermore we have D = [see 
Eq. (|A6I) ] in this limit. The critical 5k$T at which the 
splay-twist Freedericksz transition starts to prevail is ob- 
viously determined by the conditions d 2 uo(p = 0) = 
dpUo(p = 0) = at uq(p)/uf = 1- This condition is 
exploited by expanding Eq. (TT21) to order p 4 and we ar- 
rive in fact at an expression for 5k st which is identical to 
5k c in Eq. (|19p . Note that in |4| instead of exact value for 
5k c = 0.534624 the approximation 5k c 1=3 0.53 was given. 
J On the background of the exact textbook analysis pre- 
(Al Rented in this Appendix it is easy to demonstrate, why 
the authors of [T(| HH came to wrong results. Start- 
ing from the basic equations for the flexodomains in dc- 
case [equivalent to Eqs. ©] at first a set of two coupled 
ODE's that contain only second order z-derivatives has 
been obtained (see Eq. (12) in [llj]). Unfortunately the 
second step, namely diagonalizing the ODE's by a "uni- 
tary" transformation V (see Eq. (14) in 11]), does not 
work since V explicitly depends on z. Thus the main 
conclusion in [1J| that the eigenmodes associated with 
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the eigenvalues P and Q (see Eq. (20) in 11,]) would de- 
couple is untenable. This applies also to the resulting 
simple relation tan(iQd) — (see Eq. (22) in [TTJ] ) which 
would replace our correct Eq. (|12p for the neutral curve 



uq(p). Note, that the same erroneous procedure has ap- 
peared already in a previous investigation f26j | , cited also 
in [l3,[lli- 
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